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ABSTRACT 

The 3rd EGRET Catalog of High-energy Gamma-ray Sources contains 170 
unidentified sources, and there is great interest in the nature of these sources. One 
means of determining source class is the study of flux variability on time scales of days; 
pulsars are believed to be stable on these time scales while blazars are known to be 
highly variable. In addition, previous work has demonstrated that 3EG J0241-6103 and 
3EG J1837-0606 are candidates for a new gamma-ray source class. These sources near 
the Galactic plane display transient behavior but cannot be associated with any known 
blazars. Although many instances of flaring AGN have been reported, the EGRET 
database has not been systematically searched for occurrences of short-timescale (~ 1 
day) variability. These considerations have led us to conduct a systematic search for 
short-term variability in EGRET data, covering all viewing periods through proposal 
cycle 4. Six 3EG catalog sources are reported here to display variability on short time 
scales; four of them are unidentified. In addition, three non-catalog variable sources 
are discussed. 
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1. Introduction 



The Energetic Gamma-Ray Experiment Telescope (EGRET) on board the Compton 
Gamma-Ray Observatory {CGRO) began observing the high-energy gamma-ray sky in 1991. 
Since that time, EGRET has contributed enormously to the list of gamma-ray point sources; the 



third EGRET catalog ( Hartman et al. 1999 , hereafter H99) contains 271 point sources exhibiting 
significant flux above 100 MeV. Besides a single solar flare, the Large Magellanic Cloud (LMC), 
and a probable association with a radio galaxy (Gen A), the identified sources are distributed 
among two established classes of high-energy gamma-ray emitters: pulsars and radio-loud blazars. 
170 sources remain unassociated with known objects. There is great interest in the nature of 
the unidentified sources, whether they be pulsars, blazars, or representatives of new classes of 
gamma-ray sources. 

The limited spatial resolution of EGRET often hampers source identification based solely 
on position; an alternate means of determining source class is the study of fiux variability. In 
particular, pulsars are believed to be stable on short (~ 1 day) time scales ( [Ramanamurthy et al 



1997) while blazars are known to be strongly variable (e.g., Mattox et al. 1997a ). Short-term and 



long-term variability studies have also led to the discovery of two candidates for a new gamma-ray 
source class: 3EG J1837-0606 (GRO J1838-04) ( p:avani et al. 19971) and 3EG J0241-6103 (2GG 
135-FOl) (IKniffen et al. 1997| , p?avani et al. 1998| ). These sources display transient behavior but 
cannot be associated with any known radio- loud blazars. 

There are many instances of gamma-ray fiares of AGN on short time scales. For example, 
3C 279 ( IKniffen et al. 1973| , |Wehrle et al. 1998| ), PKS 0528-M34 (punter et al. 1993| ), 3G 454.3 
(IHartman et al. 1993| ), PKS 1633-^382 ( |Mattox et al. 1993| ), PKS 1406-076 ( [Wagner et al. 1995| ), 
PKS 1622-297 ( |Mattox et al. 1997a|) , and GTA 26 ( [Mattox et al. 1995| , [Hallum et al. 1997[ , [Mattox 
et al. 2000) all show strong evidence of variability on time scales of 1-3 days. Many other AGN 
show variations on time scales of weeks to months ( von Montigney et al. 1995 , [McLaughlin et aL 



1996 , Mukherjee et al. 1997 ). Although these many instances of short-term variability have been 



noted, no exhaustive survey of the EGRET database has been performed. 

It is the purpose of this study to conduct a systematic and comprehensive survey of EGRET 
data from proposal cycles 1-4 in search of evidence of short-term fiux variability above 100 
MeV. The standard method of EGRET data analysis, maximum likelihood, was applied to 
short-duration (one-day and two-day) maps of gamma-ray intensity. The variability statistic V 
was used as a guide for singling out those light curves which deviate appreciably from constant 
fiux. For the unidentified sources that display evidence of fiaring, Monte Garlo simulations were 
performed to determine the probability of such fluctuations assuming intrinsically nonvariable 
sources. Six catalog sources were found to exhibit variability on one- and two-day time scales: 
four are unidentifled, one is a fiat-spectrum radio quasars (FSRQ), and one is a BL Lac object. In 
addition, two transient non-catalog excesses and one non-catalog sub-threshold repeating source 
are discussed. The limitations of this study are then considered; a summary and concluding 
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remarks follow. 



2. EGRET Data and Preliminary Analysis 

The data used for the present analysis are those from Phase 1 through Cycle 4 of the CRGO 
mission (1991 April - 1995 October), matching the time span of the third EGRET catalog (H99). 
The CGRO timeline is divided into observation (viewing) periods at fixed attitudes that range 
in duration from one day to several weeks. Each of the 144 viewing periods longer than 3 days 
was first divided into two-day intervals (for viewing periods with an odd number of days, the last 
interval has three days). For each interval, a sky map was generated for E > 100 MeV using the 
same criteria used to generate standard EGRET maps. In addition, for each instrument pointing 
within 20° of the Galactic plane (where the sources typically are brighter), a series of one-day 
maps was generated. Each of these maps was analyzed with the EGRET maximum likelihood 
program (pVIattox et al. 1996D using the following procedure: (1) All sources from H99 within 30° 



of the pointing direction were modeled. Sources beyond 30° for these short intervals have too little 
exposure to be useful and would require separate analysis, because the EGRET Point Spread 



Function (PSF) is broader at wide angles (Esposito et al. 199£). (2) Fluxes for all the sources and 



the EGRET diffuse radiation model (Hunter et al. 1997, Sreekumar et al. 1998) within the map 



were optimized simultaneously, producing a table with statistical significance and flux for each 
source. (3) A search was made for any statistically-significant source in the map not found in the 
EGRET catalog. 

The statistics of the data are very limiting, especially for the one-day light curves. In order 
to isolate the curves for which the flux uncertainty is not overwhelming, the significance of each 
of the short-duration detections were examined. An excess of 4a is the minimum requirement for 
inclusion of a source in H99, and this threshold was applied here so that only those light curves 
which contain at least one 4(t short-duration detection were considered for further analysis. 

Even with this requirement, much of the data suffers from large statistical uncertainty. This 
leads us to the conclusion that for all but the strongest sources, EGRET is sensitive to only the 
most dramatic short-term variations in flux. In order to single out the most variable light curves, 
the variability statistic V ( McLaughlin et al. 1996| ) was used. For the remaining light curves the 



X statistic was calculated from 



■ 1 CTj- 

1=1 * 



where N is the number of observations, Fi is the detected flux during the ith observation, F 
is the mean flux for the viewing period, and ai is the 1-a flux uncertainty of the ith observation. 
If Q is the probability of obtaining a value of equal to or greater than the empirical from 
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an intrinsically nonvariable source, then V = —logQ. For this work, light curves with V > 3.0 
{Q < 10~^) are considered to manifest source variability. Light curves displaying one or more 
consecutive strong two-day detections or two or more consecutive strong one-day detections 
against a background of weak detections and/or upper limits are considered to show evidence of 
flaring. In such situations, the value of V is depressed by the high number of weak detections 
and/or upper limits; for this reason these curves were retained for further study if they had 
variability indexes at or above 1.0. 

Seven light curves of six sources displayed one or two-day variability indexes greater than 3.0 
and/or 1.0 with evidence of flaring. With one exception, the one-day analysis turned up a final 
list of variable light curves identical to that of the (independently performed) two-day analysis. 
Four of the seven sources are listed as unidentified in H99. 

For the unidentified sources that display evidence of flaring behavior, Monte Carlo simulations 
of the data were performed. The data in question were simulated using the current models of 
the diffuse isotropic radiation ( Sreekumar et al. 1998| ) and the diffuse Galactic emission ( |Huntei 



et al. 1997) in conjunction with the source flux. In some cases, the source in question was 
not isolated and a number of nearby sources were included in the simulation; all sources were 
assumed to be intrinsically nonvariable. Mean fluxes for each source and short-duration exposures 
representative of the viewing period are input into the model. For a given source, the simulation 
consisted of distributing the flux as the EGRET PSF and adding Poisson-distributed deviates 
pixel-by-pixel. In this way 1000 counts maps were generated; these maps were analyzed using 
maximum likelihood, resulting in a distribution of the significance for the source of interest. The 
probability of finding a single detection and subsequently a "flare" of given overall significance or 
greater from an intrinsically nonvariable source was then determined from this distribution. See 
Mattox et al. (1996) for more details about the Monte Carlo simulation. 



3. Discussion of Individual Catalog Sources 

The following sources are listed in H99 and display evidence for variable output. To be listed 
in II99, sources with |6| < 10° must be detected with a significance of at least 5a and sources 
with \b\ > 10° must be detected at a significance of at least Aa. Four of the sources have not 
been previously known to display short-term variability. The variability of the other two has been 
reported and for these we present additional work. 



3.1. New Reports of Short-term Variability 
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3.1.1. 3EG J0222+4253 (3C 66A/PSR J0218+4232?) 

Source 3EG J0222+4253 exhibits a strong two-day flare in VP 15.0 g; virtually all of the 
detection in this 14-day interval came on days 11 and 12. The peak flux value on day 12, 
(8.1 lb 2.4) X 10^7 photons cm~^ s~^, is over four times the mean flux for this source, either in the 
full viewing period or averaged over all observations in H99. The source was 12°5 off-axis and the 
two-day light curve (see Figure 1) registers a variability index of 2.6. Monte Carlo simulation of 
this source indicates that one should expect only about 5 in 10^ such viewing periods to display 
an equivalent or stronger "flare" if the source is intrinsically nonvariable. Source 3EG J0222+4253 
was also observed within 20° of the instrument axis in VP's 211.0, 325.0, and 427.0, but was found 
in all three cases to be nonvariable and consistent with the lower-level flux observed in VP 15.0. 

At energies above 1 GeV, the spatial analysis (H99) strongly supports an association of this 
source with the BL Lac object 3C66A {z = 0.444). This identification would place this source 
among the more than 60 blazar-class objects seen by EGRET, many of which show short-term 
time variability, as noted above. 

Verbunt et al. (1996) and Hermsen et al. (1999) suggest that much of the sub-GeV gamma-ray 
flux may have its origin in a second source, the binary millisecond pulsar PSR J0218+4232, 
located less than 1 degree from 3C66A. Their argument is based on three lines of evidence: (1) the 
100-1000 MeV data show a 3.5 sigma evidence of pulsation at the radio period; (2) the gamma-ray 
light curve resembles the one seen in hard X-rays; and (3) the spatial analysis shows that the 
source position moves from the BL Lac position toward the pulsar position with decreasing 
gamma-ray energy. 

For the two-day flare, the EGRET data do not have good enough statistics or resolution 
to determine which of the two candidate sources is responsible for the flare. Modeling both the 
AGN and the pulsar at their known positions, the likelihood analysis assigns some photons to 
each source. The flare does have many more photons at lower energies than higher energies, as 
indicated in Table 2. Also listed in Table 2 are the number of photons expected from a single 
2-day interval given the average flux of 3EG J0222+4253 in VP 15.0. From the numbers it is clear 
that the flare is seen at all energies, including the 100-300 and 300-1000 MeV bands. The energy 
distribution of the flare can also be seen in the source spectrum for all of VP 15.0, shown in Figure 
2. The photon number index is 2.37 it 0.29 for this viewing period, compared to 2.01 it 0.14 for 
the sum of phases 1-4 for this source. Because most of the photons that make up the flare are of 
lower energy, they would be associated with the pulsar and not the AGN in the two-source model. 
Flaring behavior is not seen, however, in the other EGRET-detected pulsars. We cannot rule 
out the possibility that the AGN spectrum extended to lower energies during the flare, although 
other EGRET AGN spectra tend to flatten rather than steepen during flares ( preekumar et ah 



1996i iHartman et al. 2000]) . 



^The preliminary version of this result appeared in Bloom et al. 1997. 
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3.1.2. 3EG JI4IO-6I47 

In the first four days of VP 14.0, this unidentified source decayed from a flux of almost 
(5.4 lb 1.5) X 10~^ photons cm^^ s^^ to below EGRET's sensitivity, where it remained for the rest 
of the 14-day period (see Figure 3). This is suggestive of flaring behavior. The one-day variability 
index is 1.5 and Monte Carlo simulation gives a probability of about 0.0006 that the fluctuation 
found in VP 14.0 is produced by a nonvariable source. Although the source is 27° 1 off-axis during 
VP 14.0, this was an early observation, when EGRET's sensitivity was high. 

It has been suggested that 3EG J1410-6147 (/ = 312.18, b = -0.35) is associated with SNR 
G312.4-0.4 ( ^turner, S. J. fc Dermer, C. D. 19951) , which fahs just outside the 68% error contour; 
this association is not well-confirmed. Flaring behavior would not be expected if the origin of 
the gamma rays was either cosmic rays from the SNR interacting with the surrounding material 
or a radio-quiet pulsar in the SNR. No association with any point radio source can be made 
(IMattox et al. 2000| ). Sources 3EG J0241-6103 (2CG 135-bOl; see below) and 3EG J1837-0606 



(GRO J1838-04) share this characteristic ( [Kniffen et al. 1997 , Tavani et al. 1998) of low Galactic 



latitude, variability, and lack of a strong radio counterpart. It may be that along with these 
sources, 3EG J1410-6147 represents a new class of gamma-ray emitter. This suggestion should be 
considered with caution in light of the high aspect during VP 14.0. 



3.1.3. 3EG J 1746-2851 
As this source is unidentified, strong, and coincident with the Galactic Center, it has been 



studied in some detail ( Mayer-Hasselwander et al. 1998| ). However, until now its short-term time 



history has not been examined. There is some evidence of fluctuation in VP's 16.0 and 429.0. It 
should be stated first that 3EG J1746-2851 sits in the center of the most densely-packed region of 
the gamma-ray sky; there are 10 sources listed in H99 within 10° of the Galactic Center. Given the 
broad EGRET PSF, source confusion is a serious problem. In addition, EGRET's sensitivity was 
high during this early observation but the source was 20°3 off-axis, leading to increased statistical 
uncertainty. However, while 3EG J 1746-2851 appears to fluctuate during two different viewing 
periods, no other sources in highly congested regions display any evidence of short-term variability. 

In VP 16.0 the three strongest one-day detections fall on days 7-9 of the 14-day period, 
exhibiting a peak flux of (3.40 it 0.92) x 10^^ photons cm^^ s^^ on day 8 (see Figure 4). This is 
only a 4.3(T detection, but it is flanked on days 7 and 9 by detections of 3.9c7 and 3. la respectively 
and there is no detection stronger than 2.3(T within 5 days of the peak. The one-day variability 
index for VP 16.0 is 2.1, corresponding to a probability of 0.008 that these data are consistent 
with a nonvariable source. Monte Carlo analysis is more restrictive. This source and the seven 
others within a 7° radius were modeled and there is a probability of 0.0011 that a fluctuation of 
this (or greater) significance will occur in a 14-day period given intrinsically nonvariable sources. 
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Source 3EG J1746-2851 also shows evidence of variability in VP 429.0. During this pointing 
the aspect is only 6° and the one-day variability index has a value of 3.0. The peak flux is very 
high (6.4 lb 1.7) X 10~^ photons cm~^ s~^ and on two days the source is not detected at all, but 
there is no evidence of flaring. The light curve is shown in Figure 5. 

Time variability would indicate a compact object as the gamma-ray source. Models such as 



advection-dominated accretion flow onto a massive black hole (Mahadevan et al. 1997) would be 
favored rather than a collection of pulsars or a cosmic ray enhancement. See Mayer-Hasselwander 
et al. (1998) for further discussion of possible models. 



3.1.4. 3EG J2006-2321 

During VP 13.1 (1991 October 31 - November 7), this source (/ = 18°8,6 = -26°3) exhibited 
transient behavior on time scales of 12 hours. During this VP the source was 13°6 from the 
instrument axis and the overall significance of the detection is 4.8<t. A light curve of 3EG 
J2006-2321 in VP 13.1 is shown in Figure 6. The first two days of observation are broken down 
into four 12-hour periods; the remaining points represent full days. Breakdown of the data into 
smaller time intervals is not useful as poor statistics become overwhelming. EGRET appears to 
have caught 3EG J2006-2321 on the leading edge of an intense flare, although this is not known 
with high confidence. The peak flux, centered on MJD 48560.25, is (1.75 it 0.53) x 10^^ photons 
cm~'^ s~^; this 12-hour detection has a significance of 5.4(T. The source then decays monotonically 
on this time scale; less than two days later it is not detected by EGRET at all. The flux never 
rises above 3.9 x 10"'^ photons cm~^ s~^ for the rest of VP 13.1. The ratio of peak to average flux 
for the VP is 5:1. The gamma-ray spectrum of 3EG J2006-2321 for VP 13.1 is consistent with a 
power law with exponent —2.13 it 0.31. 

Applying a test to the light curve yielded a probability of 0.0006 that these data are 
consistent with an intrinsically nonvariable source. Results from Monte Carlo simulations of the 
source are consistent with these analyses. From the simulated data it was determined that the 
probability of finding such a deviation in a 7-day viewing period is 0.0005. 

The combination of high Galactic latitude and short-term transient behavior suggests an 
association with the blazar class of AGN. While this source may well be a blazar, its radio flux 
at 5 GHz is significantly lower than blazars with similar peak gamma-ray flux values. The best 
candidate for association with 3EG J2006-2321 is the radio source PMN J2005-2310 (260 mJy 



at 5 GHz, ar not known at this frequency; (Griffith, M. R. &: Wright, A. E. 1993)). Of the 10 



high-confidence blazars listed by Mattox et al. (1997) with peak gamma-ray flux above 10^6 
photons cm~^ s~^, all are associated with flat-spectrum radio sources having a 5 GHz intensity 
above 1.0 Jy. The analysis of Mattox et al. (2000b) has determined that the a priori probability 
of EGRET detecting such a weak radio source is 0.0006. This low probability ensures that 
even though PMN J2005-2310 falls only 11.0 arcmin from the EGRET-determined position of 
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3EG J2006-2321, well within the 50% confidence countour, the probability of association is only 
0.015 ( [Mattox et al. 2000 ). However, the variability presented in the present study bolsters the 



possibility of the source being a blazar. In addition, while PMN J2005-2310 is not bright, there is 
strong evidence that it is a flat-spectrum source; it is listed with a flux of 260 mJy in the Texas 
365 MHz survey ( Pouglas, J. N. et al. 1996| ) where it is hsted as TXS 2002-233. The spectral 



index around 365 MHz is 0.7 it 0.2. The radio source also appears as NVSS J200556-231028 in the 
NRAO VLA Sky Survey ( pondon et al. 1998|) at 1.4 GHz with a strength of 302 mJy. 



A multi-wavelength investigation of PMN J2005-2310 is encouraged in order to determine 
its nature. VLBI measurements may resolve a compact core, and optical observations may aid 
in determining redshift and polarization. Finally, an X-ray detection of this source will help to 
establish its spectral energy distribution. 



3.2. Follow-up Reports of Short-term Variability 

3.2.1. 3EG JO24I+6IO3 (2CG 135+01) 

Source 3EG J024H-6103 is unidentified and has been the subject of several papers because 
of evidence of short-term variability and lack of radio-loud blazar counterparts ( [Kniffen et al. 



1997, Tavani et al. 19981). The short-term data have been re-analyzed in light of indications that 



another source (3EG J0229+6151, H99) exists only 1°68 away with about half the intensity of 
3EG J0241+6103. The results still indicate evidence of variability in VP 211.0 (1993 25 February 
- 9 March) with a two-day variability index of 3.1 (see Figure 7). Although this source is near the 
galactic plane (6 = +0°99), there are no EGRET-detected sources other than 3EG J0229+6151 
within 10°. Source 3EG J0241+6103 and the nearby source were modeled using Monte Carlo 
methods and the probability of finding a fiuctuation this strong or stronger from steady sources in 
a 14-day period is about 5.9 x 10~^. During no viewing period other than 211.0 was the source 
found to be variable. 



3. 2. 2. 3EG J 0530+1 323 (PKS 0528+1 34 ) 
Previous analyses of this quasar, located near the Galactic anticenter, found it to be strongly 



variable on two-day time scales during early EGRET observations (Hunter et al. 1993, Thompson 
et al. 1997| ). PKS 0528+134 was observed in 1993 March to have a peak flux of 3.5 x 10"'^ photons 
cm~^ s~^; by 1993 May it had decreased by a factor of 2.5. The two-day variability index of 4.6 
for VP's 0.2 through 2.1 corresponds to a probability of 0.0002 that these data are consistent with 
an intrinsically nonvariable source. This is in stark contrast to the results of almost three weeks 
of continuous observations of the Galactic Anticenter in 1995 February - March (VP's 412.0 and 
413.0). During this time PKS 0528+134 displayed virtually constant flux in its two-day light curve 
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(see Figure 8); there is no evidence of variability from this strong, previously erratic source. The 

analysis indicates that the probability that these late-period data arc consistent with a constant 
source is 0.97. Therefore one must be careful when using flux histories to determine source class; 
while strongly variable data might be used to rule out identification of pulsars, measurements of 
constant flux do not rule out identification of AGN. 

4. Discussion of Non-catalog Sources 

The thresholds for inclusion in H99 (5cr for |6| < 10° and 4a for |6| > 10°) apply to any single 
viewing period or combination of viewing periods; however, it is possible for a source to display 
a short-duration significance well above that of the entire viewing period. This is almost always 
true of flaring sources, and is true of the following two sources. Their overall detections during 
the viewing periods in question fall below the catalog thresholds, but have short-duration excesses 
above these values. 

4.1. GRO J0927-41 

This source appears near I = 266.68, b = 7.03 with a 95% confidence radius of 1°01 during 
the third and fourth days of VP 338.5. On these days the source displays excesses of 3.3(7 and 
4.0(7 respectively (see Figure 11). The combined significance for these two days is 5.2(7. The very 
strong Vela pulsar is 10°4 away; Monte Carlo simulation of GRO J0927-41 and Vela indicates 
that the probability of such transient behavior from an intrinsically nonvariable source is 0.0013. 
A search for counterparts revealed PKS 0920-397, a FSRQ 1°33 away near the 99% confidence 
contour. This quasar (z = 0.591) has a flat radio spectrum above 1.4 GHz, with a flux of 1.5 Jy at 
4.85 GHz. It may be that EGRET has detected the tip of a flare from this source. There are not 
enough counts in the data to produce a spectrum of this source. 

4.2. GRO J1547-39 

In VP 226.0, an excess of 4.7(7 appears near I = 336.17, b = 11.78 95% with a confidence 
radius of 0°97. The source is given the name GRO J 1547-39 and is seen most clearly in the third 
of four two-day maps. The light curve is shown in Figure 10. When broken down into one-day 
intervals, the significance falls to 4. Oct and 3.1ct. The overall detection for this source in VP 226.0 is 
3.6(7; in no other viewing periods that include the source position does the overall significance rise 
above 3.0(7. No flat-spectrum radio loud blazars or other candidate objects appeared as candidates 
for association with this source, and there are not enough counts to provide a spectrum. Monte 
Carlo simulation of this source indicates that there is a probability of 0.010 of finding a "flare" of 
equal or greater strength from an intrinsically nonvariable source. With such a high probability 
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we make only a weak claim for variability of GRO J1547-39. 



4.3. Sub-threshold Repeating Sources 

In the region —40 < b < 40, 148 two-day non-catalog excesses with significance > 3.5cr are 
found. Among these are seven spatially coincident pairs and one spatially coincident triplet. 
For each set, the data are combined; in six of these combined maps the significance exceeds the 
standard thresholds. These pairs are considered for further analysis, with the hypothesis that each 
one represents the same transient source at high flux values. 

This hypothesis is almost certainly not true. The probability of finding six such coincident 
pairs from 148 excesses is not negligible; assuming that the sources are randomly scattered and 
that each source has an error circle with radius ~1°0, the probability of finding six or more 
spatially coincident pairs is nearly 0.03. Beyond this simple analysis, the individual excesses are 
below threshold and most are close to the Galactic plane, and may not represent real sources. 
However, of these six pairs one is worthy of mention for its high combined significance. 

During VP 423.5, a source appeared near I = 334.78, b = —0.28, prompting a report to the 
lAU ( IKanbach et al. 1995| ). The source, GRO J1627-49, has a 95% confidence radius of 0°42 and 



was found to be coincident with an H II region. However, it was not seen in any other band and 
faded during the viewing period; it was not strong enough for the viewing period as a whole to be 
included in H99. The present analyses reveal that during VP 23.0 a sub-threshold excess appeared 
for two days at a coincident position. The data from these viewing periods combine to yield an 
excess with significance 5.7a. The spectrum of these combined data is shown in Figure 11 and is 
rather flat, with a spectral index of 1.78 it 0.31. 



5. Limitations of this Study 

egret's sensitivity, while higher than any previous gamma-ray telescope, is low enough 
to severely limit the present analyses when combined with the intrinsically low gamma-ray flux. 
Often there are simply not enough counts to break t he data down into smaller divisions and retain 
any meaningful information. For much of the data, low counts translate into statistical errors of 
~40-50%. In these cases, an intrinsic flux increase of 100% is barely detected with any significance; 
changes of ~200-300% are required for higher-confidence claims of variability. Thus only the most 
dramatic changes in flux are measured with any confldence. For all but the strongest sources, the 



estimated 10% systematic uncertainty (Thompson et al. 1995) adds insigniflcantly to statistical 
error. 



Interpretations of the variability index V are rendered problematic by the low EGRET count 
rate. A low value of V may have its origin in constant source output; it may also be due to poor 
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statistics. There is no way to distinguish between these two cases. Again, this means that only 
very large values of V are meaningful, and that only strong variations from constant flux are 
detected. 

Also, the high short-duration errors lead to uncertainties in mean fluxes for given viewing 
periods; the shorter the duration of the detection, the higher the uncertainty of the mean. This 
may result not only in misleading values of V but also Monte Carlo probabilities, as these 
calculations assume the averages for the viewing period in question. 

6. Summary and Conclusions 

The results of this work are summarized in Table 1. The first column indicates the source 
name by increasing Right Ascension; sources found in H99 are given the 3EG designation and 
non-catalog sources are given the GRO designation. The second column indicates alternate names, 
and the third column indicates source identification as FSRQ, BL Lac, or unidentified (UID). 
Column four indicates those viewing periods in which the source is found to be variable, column 
five gives the variability index for that viewing period, and column six indicates whether the 
source flared in that viewing period. The Monte Carlo probability for unidentifled flaring sources 
is given in the seventh column. 

The detections of short flares in 3EG J0222+4253 and 3EG J2006-2321 demonstrate that 
EGRET can detect short, intense flares in relatively weak catalog sources (3.9(7 and 4.8a in the 
viewing periods in which the flares were seen) as well as in strong blazars. Therefore most EGRET 
sources, except bright AGN, are not characterized by strong variability on one- and two-day time 
scales. If intense flares from sources like 3EG J0222-I-4253 and 3EG J2006-2321 were common, 
this survey would have detected more of them. 

The same conclusion is reached for sub-threshold sources. The detections of GRO J0927-41 
and GRO J1547-39, plus the repeating sources, demonstrate that EGRET can find very weak 
transients from sub-threshold sources with this sort of systematic analysis. The fact that few such 
intense flares are seen in this study indicates that the gamma-ray sky is not flUed with these kind 
of flares from sub-threshold sources. 

The flaring behavior seen in 3EG J0222-I-4253 adds another puzzle to the issue of identiflcation 
of this source. The flaring is seen at energies between 100 MeV and 1 GeV, where the emission 
may be dominated by PSR J0218-I-4232. This would make this source unusual in two ways; it 
would be the only ms pulsar seen by EGRET and the only pulsar to show strong flaring. 

The long observation of steady emission from PKS 0528-1-134 shows that, as seen by EGRET, 
blazars are not always variable even on time scales of three weeks. Therefore long episodes, free of 
intense flares from a given source, do not rule out the identification of that source as a blazar. 

The observation of clear short-term variability in 3EG J2006-2321 indicates that this source 
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may be a blazar, although its radio flux at 5 GHz is significantly less than other EGRET-detected 
blazars with similar peak gamma-ray brightness. The suggestion of variability of 3EG J1410-6147 
may place it alongside 3EG J0241-6103 and 3EG J1837-0606 as a third candidate for a new 
gamma-ray class of variable Galactic sources. 

In summary, the results of this study suggest that EGRET reached the threshold of an 
important new feature of gamma-ray astrophysics: the ability to measure short-term time 
variability. The few examples found by EGRET show that such sources do exist and axe not 
restricted to the blazar class. Variability offers an effective tool for matching gamma-ray sources 
to objects seen at other wavelengths. With its ability to detect much fainter sources, measure 
variability on time scales of hours, and detect smaller flux variations that EGRET, the proposed 
Gamma-ray Large Area Space Telescope (GLAST) should be able to take full advantage of 
gamma-ray source time variability to extend the results of the present study. 

The authors wish to thank J. R. Mattox and J. P. Halpern for their helpful comments and 
suggestions. P. W. gratefully acknowledges support from the NASA/ASEE Summer Faculty 
Fellowship Program. 
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Fig. 1.— Light curve of 3EG J0222+4253 (3C 66A) during VP 15.0. The points represent 48-hour 
integration times. The variabihty index V = 2.6. 

For all light curves, the dashed line indicates the average flux for the viewing period. 
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Fig. 2.— Spectrum of 3EG J0222+4253 (3C 66A) during VP 15.0. The spectral index is 2.37±0.29. 
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Fig. 3. — Light curve of 3EG J1410-6147 during VP 14.0. The points represent 48-hour integration 
times. The variabihty index V = 2.4. 
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Fig. 4. — Light curve of 3EG J1746-2851 during VP 16.0. The points represent 24-hour integration 
times. The variabihty index V = 2.1. 
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Fig. 5. — Light curve of 3EG J1746-2851 during VP 429.0. The points represent 24-hour integration 
times. The variabihty index V = 3.0. 
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Fig. 6. — Light curve of 3EG J2006-2321 during VP 13.1. The first four points represent 12-hour 
integration times and the last five points represent 24-hour integration times. The dashed hne 
indicates the average flux for the viewing period. The variabihty index V = 3.2. 
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Fig. 7.— Light curve of 3EG J0241+6103 (2CG 135+01) during VP 211.0. The points represent 
48-hour integration times. The variability index V = 3.1. 
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Fig. 8.— Light curve of 3EG J0530+1323 (PKS 0528+134) during VP's 412.0 and 413.0. The 
points represent 48-hour integration times. The variability index V = 0.013. 



-23- 




49595 49600 49605 49610 

JD - 2400000.5 



49615 



Fig. 9. — Light curve of GRO J0927-41 during VP 338.5. The points represent 24-hour integration 
times. 
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Fig. 10. — Light curve of GRO J1547+42 during VP 226.0. The points represent 48-hour integration 
times. 
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Fig. 11.— Combined spectrum of GRO J1627-49 from VP's 23.0 and 423.5. The spectral index is 
1.78 ±0.31. 
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Table 1. Sources Reported in this Study. 



Source 


Alt Name 


ID 


VP's 


V 


Single Flare 


p 


3EG J0222+4253 


3C 66A/ 


BL Lac/ 


15.0 


2.6 


yes 


0.0005 




PSR J02 18+4232 


PSR 










3EG J0241+6103 


2CG 135+01 


UID 


211.0 


3.1 


yes 


0.0006 


3EG J0530+1323 


PKS 0528+134 


FSRQ 


0.2-2.1 


4.6 


no 


n/a 


GRO J0624-02^ 


none 


UID 


1.0,413.0 


n/a 


yes 


n/a 


GRO J0927-41 


none 


UID 


338.5 


0.5 


n/a 


0.001 


3EG J1410-6147 


none 


UID 


14.0 


2.4 


yes 


0.0006 


GRO J1547-39 


none 


UID 


226.0 


1.4 


yes 


0.010 


GRO J 1627-49'' 


none 


UID 


23.0,423.5 


n/a 


n/a 


n/a 


3EG J1746-2851 


none 


UID 


16.0 


2.2 


yes 


0.001 








429.0 


3.0 


no 


n/a 


3EG J2006-2321 


none 


UID 


13.1 


3.2 


yes 


0.0002 



Sub-threshold repeating source. Light curves were not generated for this source since it did 
not display any short-duration excess greater than 4cr. 



Tabic 2. Fhixcs for the Two-day Flare of 3EG J0222+4253. 



Energy Range Significance Flux Photons Expected Photons 



100-300 MeV 3.2a 35 ± 14 18 ± 7 6 ± 2 

300-1000 MeV 4Aa 22 ±8 13 ±5 3±1 

> 1000 MeV 1.8a 4 ±3 1.5 ±1.4 0.2 ± 0.2 



